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Aquarius Salinity Retrieval Algorithm: Final Pre-Launch Version 

Frank J. Wentz and David M. Le Vine 


1 Introduction 

This document provides the theoretical basis for the Aquarius salinity retrieval algorithm. 
The inputs to the algorithm are the Aquarius antenna temperature (Ta) measurements along with 
a number of NCEP operational products and pre-computed tables of space radiation coming from 
the galaxy and sun. The output is sea-surface salinity and many intermediate variables required 
for the salinity calculation. This revision of the ATBD is intended to be the final pre-launch 
version. 

The ATBD is organized into two main sections. The Forward Model Section describes how 
the Aquarius antenna temperature is calculated given the following information: 

1. Footprint location: time, latitude, and longitude 

2. Pointing angles: Earth incidence and azimuth angles; sun and moon pointing vectors 

3. Antenna pattern 

4. NCEP operational fields: surface temperature and vector wind, atmospheric profiles 

(pressure, temperature, humidity, liquid water), sea ice 

5. Solar flux values from radio astronomy observations 

6. Pre-computed tables giving the Ta contribution of the galaxy and sun (moon contribution 
is computed from an analytic expression) 

The Retrieval Algorithm Section describes how the various components of the forward 
model are used to estimate salinity given the Aquarius Ta measurement. The algorithm is 
essentially a subtraction process in which the unwanted sources of radiation (galaxy, sun, moon, 
and Earth’s atmosphere) are removed from the Ta measurement in order to obtain just the 
Earth’s surface emission term. A simple regression is then use to estimate salinity from the 
surface emission. 


2 The Forward Model 

The forward model gives the Aquarius antenna temperature in terms of the antenna pattern 
G convolved with the surrounding brightness temperatures (Tb) 47i-steradian field. The 
convolution of G with Tb is called the antenna temperature equation. The field of surrounding 
brightness temperatures is partitioned into a space segment and an Earth segment. 


2.1 The Antenna Temperature Equation 

In this subsection, we work in terms of the classical, rather than the modified, Stokes 
parameters. The measurement vector Ta , mea , which is called the antenna temperature vector, is 
then defined as: 
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where 7 \mea,Vt T Amea.Hi T Amea. +45, T Amea. -45 - T Amea.left, and T Amea, right ^Te the V pol, h pol, plus 45 , 

minus 45°, left circular, and right circular polarization measurements. Aquarius will only 
measure T Ame a,v, T Am ea,H, T Am ea,+45 > and T Amea , -45 1 but the full Stokes formulation is included here 
for completeness. Likewise, brightness temperatures are defined as 
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The measurements represent an integration over the entire 4k steradians surrounding the 
antenna. We divide this integration into 2 components: the Earth field of view and the space 
field of view. 


T = T +T 

A ,mea A, earth A, space 

(3) 


(4) 

Earth ^ 


space = | G (b)T I{ s/wrp (7Q 

(5) 


Space 


The first integral is over the surface of the earth, where the differential surface area is clA. The 
second integral is over space, where <r/Q is the differential solid angle. The other terms in these 
equations are as follows. The matrix G is a 4x4 matrix of the antenna gain function. Each 
element in this matrix is a function of the look direction b. For the first integral, b is the unit 
vector pointing from the antenna to dA. For the second integral, b is the unit vector in the 
direction specified by dil. 


The term 'P(^) is a rotation matrix defined by 


I'M 


10 0 0 
0 cos 2 (f> -sin 2(f) 0 
0 sin 2(f) cos 2<f> 0 

0 0 0 1 


(6) 


The rotation angle </> is the sum of the antenna polarization angle <f> ant and Faraday rotation angle 
(j)f due to the ionosphere: (j)=(j)ant+(j)f . The angle <j) an , is the angle between the antenna polarization 
vectors and the Earth polarization vectors. It is an implicit function of b as well as the attitude 
of the Aquarius spacecraft. The Faraday rotation angle is a function of b as well as the 
Earth’s magnetic field vector B and the electron density along the path b. 
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where v is frequency in GHz, N e is vertically-integrated electron counts in TEC units, B is the 
Earth’s magnetic field vector in nanotesla units at the mean height of the ionosphere, and the last 
term is the partial derivative of slant range to the vertical height, which converts N e from a 
vertically-integrated value to a slant-range integrated value. The rotation matrix is not required 
for the space integration because we assume the space radiation is unpolarized. 

The ratio of the differential solid angle to the differential surface area is 


dCl_ cos#, 

^ A J lat 2 

oA r 


( 8 ) 


where 0, is the incidence angle and r is the range. For a spherical Earth, the leading term ()„, 
would be unity. However, the Earth is modeled as an oblate spheroid and as a consequence this 
term is a function of latitude, deviating about +1% from unity. 

The brightness temperature vector for Earth in equation (4) is 
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where T B v,toa and Tku mm are the v-pol and h-pol Earth brightness temperatures measured at the 
top of the atmosphere (TOA), but below the ionosphere. We have assumed that the 3 ld and 4 th 
Stokes parameters coming from the Earth are negligibly small at 1.4 GHz. 

The brightness temperature vector for space is 
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2 T, 
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( 10 ) 


The space contribution consists of cosmic background radiation, galactic radiation, and radiation 
from the sun and the moon. We assume the space radiation is unpolarized (i.e., only the 1 st 
Stokes is non-zero). (There may be a small 4 th Stokes component, which needs to be 
investigated.). Section 2.2 discusses this space contribution in more detail. 

The Earth’s limb also contributes to the radiation received by the antenna. The Earth’s limb 
varies from 260 K at the surface to 11 K at 10 km. Above 40 km, the limb brightness 
temperature is less than 0.001 K. When integrated over the antenna pattern, the limb contributes 
only 0.006 K to the 1 st Stokes and is essentially zero for the 2 nd through 4 th Stokes. 

At this point in our analysis, we drop the 4 th Stokes from the notation. Aquarius does not 
measure the 4 th Stokes, and it does not need to be considered. So hereafter, all T A and T B vectors 
now have just 3 components, and the antenna gain and rotation matrix are 3 by 3. 
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2.2 Radiation from Space 

Space radiation consists of cosmic background (2.73 K), galactic, solar, and lunar 
components. This radiation is received by the antenna in two ways: directly and via Earth 
reflection and scattering. The contribution of direct lunar radiation is completely negligible 
(lunar radiation is a factor of 10 4 less than the solar radiation). However, lunar radiation 
reflecting off the ocean surface is not negligible because at certain times each month the 
reflected ray enters the antenna mainbeam. 

These space radiation terms are denoted by: 

1. Direct and reflected galactic radiation: Ta , ga i_dir and Ta , S ai_ref 

2. Direct, reflected, and backscatter solar radiation: Ta, sun_dir ? t a , sun_ref > and Ta, sunjbak 

3. Reflected lunar radiation: Ta. ,mon_ref 

Because the solar radiation is so intense, one must consider both the specular reflected 
component that enters the far sidelobes of the antenna and the backscatter component that enters 
the mainbeam of the antenna. 


Given these 6 terms, equations (4) and (5) can be partitioned as 

T = T + T + T + T + T 

A, earth A, earth _dir A, gal _ ref A, sun _ ref A,sun_bak A,mon_ref 


( 11 ) 


A, space A,gal_dir A,sun_dir ( ) 

where the new term Ta , ea rth_dir refers to radiation coming just from the earth (i.e., no space 
radiation reflections). It is better to directly compute the of antenna temperatures of the space 
reflections by numerically performing the required integrations rather than trying to model them 
as brightness temperature components reflected off the earth. For example, one could model the 
reflected galactic radiation as one component of T\\ Jna in (4). However, the antenna pattern 
correction discussed in Section 3.3 works best for a smoothly varying brightness temperature 
scene. The galactic radiation field has very sharp features that would present a problem for the 
APC. Therefore, we decided to numerically compute the antenna temperatures for the various 
space reflections and store the results of the integrations as tables. The next subsections describe 
how the 6 space terms are computed. 


2. 2 . 1 Direct Galactic Radiation 

The direct galactic radiation Ta ga i_dir is computed by numerically evaluating the integral in 

(5), 


T*^.*= £ | G (b)T„. w ,dn 

Space 


(13) 
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2T ; 
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0 

0 


(14) 


where T B ga i comes from a galactic map [Le Vine and Abraham, 2004] that has the cosmic 
background radiation of 2.73K added in. The pixel size in the LeVine-Abraham map is 
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0.25°x0.25°, and the minimum pixel value is about 3K. We interpret the 3 K floor of the galactic 
map as the cosmic background radiation (2.73 K) plus radiation coming from distant galaxies. 

The computation of Ta , ga i_dir depends on the orbit position of the satellite and on the time of 
year as the Aquarius orbit precesses relative to the fixed galaxy every sidereal year (365.25636 
days). A separate computation is done for each of the 3 radiometers because the antenna patterns 
are different. Tables of Ta , ga i_dir are generated for 1441 positions within an orbit (every l A°) and 
for 1441 periods during a sidereal year (about every 6 hours). These tables are specified by 
tagal_dir_tab(1441 ,1441 ,3,3), where the 4 dimensions correspond to time of sidereal year, orbit 
position, Stokes number, and radiometer number. Time is referenced to January 1, 2010, 0Z, and 
orbit position is referenced to the South Pole node. The South Pole node is the point in the orbit 
where the polar component of spacecraft velocity vector changes from pointing south to pointing 
north. The orbit position is the angle between a vector pointing from the Earth center to the 
South Pole node and a vector pointing from the Earth’s center to the current position of the 
satellite. For example, 90° would correspond approximately to the satellite crossing the equator 
going south to north. The periodicity of the time variable is one sidereal year (365.25636 days). 
Even though we are assuming the galactic Tb is unpolarized, antenna polarization mixing results 
in Ta , ga i_dir having all four Stokes components, although we drop the 4 th Stokes. Ta , ga i_dir is a 
small (0.3K), slowly varying term and the granularity size of the table is more than sufficient. 
During operational processing a bi-linear interpolation is used to look up tagal_dir_tab as a 
function of orbit position and time of year. 

The calculation of Ta , ga i_dir depends upon the assumed ascending equator crossing time for 
Aquarius. A value of 5:59:59.16 pm was used to generate the table tagal_dir_tab. If the actual 
ascending node time is different, the relative phase of the table can be adjusted via an input 
constant orbit_phase_dif that specifies the actual node time relative to 6:00:00 pm. Currently 
orbit_phase__dif is set to -0.84 seconds, and it will be revised after launch to the true value. 


2.2.2 Reflected Galactic Radiation 

The reflected galactic radiation TA, ga i_ref is the most difficult space term to deal with. It can 
be large (5 K), have sharp features, and be modulated by the ocean surface roughness. To 
calculate Ta ,gai_ re f, an integration over the Earth is done as follows 

J G(b)fWV^ (15) 

Earth 

where Tb , ga i_ref is the brightness temperature of the reflected galactic radiation at the top of the 
atmosphere and is given by 

^15 ,gal_ref =T ( R Bgal ~ R Bcos ) R (^) 


R = 


R v + R h 

R v - R h 

0 


(17) 


where z is the atmospheric transmittance and vector R is the reflectivity for the 1 st and 2 nd 
Stokes in terms of the v-pol and h-pol reflectivity. Equation (16) shows that we subtract the 
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constant Tbcos = 3K from the LeVine -Abraham map before computing the reflected galaxy 
radiation. This is done because T Bc os is the one and only space radiation term that is included in 
Earth radiation calculation (see Section 2.3). Hence, we must subtract it here to avoid including 
it twice. 


Note that (16) only accounts for galactic reflections in the specular direction. In truth, 
bistatic scattering from a rough ocean will scatter galactic radiation from many different 
directions into the mainlobe of the antenna. In effect, a rough ocean surface tends to add 
additional spatial smoothing to Ta , g ai_ref- By neglecting this scattering effect we introduce a 
significant deficiency in the current forward model and retrieval algorithm. It is a complicated 
process that will be addressed between now and launch. The plan is to replace the current 
galactic reflection tables with new tables that include the scattering effect. 

Tables of Ta , ga \_ref are made. As was the case for the direct galactic radiation, Ta , g ai_ref 
depends on the orbit position of the satellite and on the time of year. The tables have the same 
structure and format as the Ta , ga i_dir except that one additional dimension is included for wind 
speed. The tables have the form tagal_ref_tab(1 441 ,1 441 ,3,3,5), where the first four dimensions 
are the same as described above for Ta , ga i_dir- The 5 th dimension is wind speed going from 0 to 
20 m/s in 5 m/s steps. 

When computing tagal_ref_tab(1 441 , 1 441 ,3,3,5), the atmospheric transmittance z is set to 
unity and the reflectivities Ry and Ru are set to ocean reflectivity values corresponding to a 
surface temperature and salinity of 20°C and 35 psu, respectively, and wind speeds equaling 0, 5, 
10, 15, and 20 m/s (i.e., the last dimension of the table). Also, the Faraday rotation angle is set 
to 0 when computing (15). During operational processing, an adjustment to the table values is 
made to convert these nominal values to the actual values corresponding to a given observation. 
This adjustment is now discussed. 

During operational processing, NCEP and other ancillary data are used to determine the 
actual values for z and R v and R H . Also, the Faraday rotation angle (j)f is computed as discussed 
in Section 3.4. The T^ ga i_ re f table values are adjusted to correspond to these actual values z, Rv, 
Rh and (f>f as follows. First, the T\. ga i_r e f table values are converted to brightness temperatures 
using the APC equation (46) below. These classical Stokes Tbs are converted to conventional v- 
pol and h-pol Tb. Then the following adjustment is made: 


T' 

1 BP,gal_ref 



Rp_ 

R 0P 


T 

± BP,gal_ref 


(18) 


where the prime sign denotes the adjusted value and subscript p denotes polarization ( V or H). 
The value r is the NCEP transmittance, Rop is the reflectivity computed at the nominal of 20°C 
and 35 psu, and Rp is the reflectivity computed at the NCEP temperature and reference salinity. 
T bp g ai ref > s then converted back to classical Stokes, Faraday rotation is applied, and the inverse 

APC equation (i.e. A" 1 ) is applied to convert back to antenna temperature. Simulations showed 
that there was less error in doing the adjustment at the Tb level as compared to simply applying 
the adjustment to the TaS. 
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2. 2. 3 Direct Solar Radiation 

The direct solar radiation Ta , sun _dir is computed in the same way as for Ta , ga i_dir- However, 
the evaluation of the integral in (5) is much simpler because the sun is a very localized source 
and can be removed from the integral. Doing this gives 


T =T 

A,sun_dir B,sun 


Q. 


sun / U \ 

^ ^ w sun _ direct ) 


2 

0 

0 

0 


(19) 


where Q su „ is the solid angle subtended by the sun and b sun _direct is the unit vector pointing from 
the spacecraft directly to the sun. The sun’s brightness temperature is given by 


T 

B.sun 


a 2 F 
2 kn 


( 20 ) 


where X is the radiation wavelength, F is the solar flux, k is Boltzmann’s constant, and f2 sun is 
about 8.216E-5 steradians. This assumes an angular radius of 0.293°, which is 10% greater than 
the optical angular radius [Wentz, 1978]. 

As with the direct galactic radiation, the computation of Ta , sun _dir depends on the orbit 
position of the satellite and on the time of year. Tables of Ta , sun _db- are generated in the same 
way as for Ta , ga i_dir and have the same structure and format, namely tasundir_tab(1441 ,1441 ,3,3). 
That is to say, tagal__dir_tab(1 441 ,1441 ,3,3) and tasun_dir_tab(1441 ,1441 ,3,3) are completely 
analogous with one exception. For the sun table, a change in the ascending node time cannot be 
simply handled by an orbit phase change. Rather the table will need to be regenerated after 
launch (a quick process that takes about an hour) when the true ascending node time is known. 

In generating tasun_dir_tab(1441 ,1441 ,3,3), we set the flux value F to 1 solar flux. Then for 
operational processing, the table values are multiplied by the actual flux value obtained from 
radio astronomy measurements. The solar flux F is a standard component of Aquarius ancillary 
data. 

An analysis of the effect of solar radiation on Aquarius is given in Wentz [2005, 2007], 
which provides maps showing how the solar contamination varies during the course of a year. 
Typical values for Ta , sun _dir are 0.05 K or less. However, because our knowledge error for 
antenna gain G in the sun’s direction may be large, the calculation of Ta , sun _dir is subject to error. 
Post-launch analyses will be done to determine if Ta , sun _dir is indeed small enough to neglect. 
The processing algorithm contains an option for setting Ta , sun _dir to zero. 


2.2.4 Reflected Solar Radiation 

The reflected solar radiation T\. v „„_,. f ,/- comes from the specular reflection of sunlight from a 
location on the Earth far away from the observation cell. The reflected sunlight enters into the 
far sidelodes of the antenna. Even though the gain of these far sidelobes is very small, the 
intensity of the sun at 1.4 GHz is so large as to make this a potential problem. The reflected 
solar component can be found explicitly from the above formulation to be (also see Wentz [1978] 
and LeVine et al. [2005], Dinnat and Le Vine [2008]) 
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<21) 

where the prime sign on the sun solid angle Q sun indicates that the apparent solid angle of the sun 
reflecting off the spherical Earth is less than the solid angle when viewed directly. This effect is 
analogous to viewing an object through a convex mirror. The object will look smaller than it 
really is. This effect is discussed more in Section 2.2.6. The brackets in (21) denote an average 
over the neighborhood centered on the location of specular reflection, which is far removed from 
the observation footprint. The amount of averaging depends on surface roughness. For a 
specular (i.e., perfectly flat) surface, no averaging is required, and (21) can be evaluated for b 
pointing in the direction of the specular reflection. For high winds when the surface becomes 
rough, the antenna gain needs to be averaged over 5° to 10°. In the sidelobe region, the terms 
T 2 , Th (j>), and R vary more slowly than G(b), and (21) can be approximated by 

=T,„ ^“-< f G(b)'r(**)R,* (22) 

where the subscript ref denotes the quantity is calculated assuming a perfectly specular reflection 
and the overbar on the antenna gain function denotes that the antenna gain has been smoothed 
over a angular range to account for surface roughness. 

The calculation of Rv.ref, Rti.ref, and T re f at the remote sun reflection point adds additional 
complexity to the modeling and algorithm. All of the ancillary data are keyed to the observation 
footprint which can be thousands of kilometers away from the sun reflection point. Furthermore, 
simulations [Wentz, 2005, 2007] indicate T\ Mm _ re f is very small, 0.01 K or less, and can probably 
be neglected. For now, we simply assume the reflectivity is that for a specular ocean surface at a 
temperature of 20°C and a salinity of 35 psu and the transmittance is set to unity. If post-launch 
analyses indicate Ta , sun _ref has a measureable effect on the observations, then the antenna gain 
function G(b) will need to be revised (i.e., increased in the sidelobe region) and the retrieval 
algorithm will need to be modified to compute the reflectivity and atmospheric transmission at 
the remote location where the sun’s specular reflection occurs. 

Tables of Ta.™, ,_,•<>/ are generated in the same way as for Ta , SU n_dir and have the same structure 
and format, namely tasun_ref_tab(1 441 ,1 441 ,3,3). The processing algorithm contains an option for 
setting T\' S un_refl to zero. 


2.2.5 Backscattered Solar Radiation 

Most of the time, the Aquarius footprints will be entirely in the dark. However at the higher 
latitudes in the summer, the sun will be visible near the horizon, and solar radiation will impinge 
onto the Aquarius footprint at a grazing incidence angle. A very small portion of this sunlight 
will be backscattered into the mainlobe of the antenna. We denote this backscatter component as 
Ta ,sun_bak, and it has a magnitude on the order of 0.1 to 0.2 K when it does occur. 

We use a model developed by Dinnat and Le Vine [2008] to specify Ta , sun _bak- This model 
parameterizes Ta , sun _bak in terms of the sun’s zenith angle and wind speed. Ta , sun _bak is 
tabularized as tasun_bak_tab(161 ,26,3,3), where the 4 dimensions corresponds to sun zenith angle, 
wind speed, Stokes number, and radiometer number. The zenith angle goes from 58° to 90° in 
steps of 0.2°, and the wind speed goes from 0 to 25 m/s in steps of 1 m/s. This table was 
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generated assuming a solar flux F of 264 solar flux units. For operational processing, the values 
from tasun_bak_tab(1 61 ,26,3,3) are multiplied by F/264, where F is the actual flux value obtained 
from radio astronomy measurements. 

The Ta ,sun_bak tables are computed assuming a nominal surface temperature and salinity of 
20°C and 35 psu, an atmospheric transmittance of unit, and no Faraday rotation, i.e., the same 
nominal conditions used for the TA, ga i_ref tables. During operational processing, an adjustment to 
the table values is made to convert these nominal values to the actual values corresponding to a 
given observation. This adjustment is exactly the same as described above for the reflected 
galactic radiation, i.e., equation (18). It should be pointed out that since Ta , sun _bak is a very small 
term (0.1 to 0.2 K), the adjustment done by (18) has negligible effect over the ocean. The only 
appreciable effect of the adjustment is that over land the ratio Rp/Rop becomes very small (land 
has a much smaller reflectivity than ocean), thereby reducing the value of Ta , sun _bak, which is the 
correct thing to do. 


2. 2 . 6 Reflected Lunar Radiation 

At certain times each month, moonlight reflected off the ocean surface enters the mainlobe 
of the antenna. This phenomenon is discussed by Dinnat et al. [2009] and is modeled here by 
using a simplified antenna gain function. A simplified approach can be used because the effect 
occurs rarely and is small in magnitude (1 K or less). The simplified gain pattern for the 1 st and 
2 nd Stokes is 


G = 


G n 

V G 21 




7 22 J 



(23) 


where <f is the angle between the antenna boresight vector and the vector from spacecraft to 
moon specular reflection point. The angle do is the half-power angle of the antenna pattern. Its 
values for the inner, middle, and outer horns are 3.04°, 3.17°, and 3.24°. The values for the 
boresight gains Gij are given in the computer code for each horn. Equation (23) is a very 
commonly used approximation for the gain over the mainlobe of the antenna. Given the small 
size of the effect, this approximation is perfectly adequate. 

The reflected lunar radiation is then given by 


T 

A ,mon _ ref 


1 B,moon 


J=-r 2 GR 
4 n 


(24) 


where the reflection vector R is given by (17) but in this case the last component, which is zero, 
is discarded in order to match the 2x2 rank of G . The terms R and r are computed from the 
ancillary data for the observation. T B , m0 on is assumed to be 275K and Q' m00 n is the apparent solid 
angle of the moon’s reflection off the spherical earth. The relationship between the true solid 
angle Q when viewed directly and the apparent solid angle Q' when viewed as a reflection is 

° ,2 ( 25 ) 

I 2s i 
1+— 

V V 
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where s is the slant path from the satellite to the point of reflection, and Re is the Earth’s radius. 
£V for the inner, middle, and outer horns is 3.93e-05, 3.79e-05, 3.63e-05 steradians. No 
adjustment is made for Faraday rotation, which is a small effect for the reflected lunar radiation. 


2.3 Earth-Only T B Equation 

The one remaining term to be specified in the forward model is T \^ a nh_dir, which is the 
radiation coming just from the earth (i.e., no space radiation either direct or reflected). The one 
exception, the uniform cosmic background radiation that reflects off the Earth, is included in 

Ta ,earth_dir 

j G(b)'P( ( i)T M! ^d4 (26) 

Earth 

where it is understood that Tbe refers to the TOA brightness temperature for radiation coming 
just from the earth. When doing forward modeling simulations, the integral in (26) is found by 
doing a very precise integration over the Earth’s surface at a spatial resolution of about 1 km 
over the mainlobe of the antenna (a coarser resolution is used outside the mainlobe). This 
section discusses the Earth-only brightness temperature Tbe that appears in the kernel of the 
integral. 


2.3.1 Theoretical Formulation 

The AMSR-E ATBD [Wentz and Meissner, 2000] gives a very detailed description of the 
radiative transfer equation for Tbe- This formulation is also applicable to Aquarius. In this 
section, it is more convenient to work in terms of the modified Stokes parameters (i.e., v-pol and 
h-pol) rather than the classical Stokes parameters. We are assuming that at 1.4 GHz, T B e 3 ld and 
4 th Stokes are zero. The v-pol and h-pol components of T B e are given by 

T be.p =T bu + t [EpTs + T Bpn ] (27) 

where subscript P denotes polarization, either V or H. Tbu is the brightness temperature of the 
upwelling atmospheric radiation, r is the atmospheric transmittance, Ep is the Earth surface 
emissivity, T s is the Earth surface temperature (K), and T B pn is the downwelling sky radiation 
that is scattered off the Earth surface in the direction of the observation. 


The atmospheric terms can be expressed in terms of the transmittance function z(si,S 2 ) 

f A 


r (-SpSj) = exp - J ds a(s) 


v * 


(28) 


which is the transmittance between points sj and S 2 along the propagation path 5. The term a(s) 
is the atmosphere absorption due to oxygen, water vapor, and liquid water. See Wentz and 
Meissner [2000] for a more complete description of a(s). The total transmittance r in (27) is 
given by 

r = r(0,S) (29) 
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where S is the slant range from the surface to the top of the atmosphere. The upwelling and 
downwelling atmosphere emissions are given by 

T bu = J<A a(s ) T(s ) t(s,S ) (30) 

0 

5 

T BD = ^dsa(s) T(s ) r(0,s) (31) 

o 

where T(s ) is the air temperature profile. The scattering integral for the downwelling radiation is 

7u/ 2 2n 

T Bpn j sin^J^(r BZ) +rr B cos)[o- 0iC (k s ,ki) + o- () x(k s ,k i )] (32) 

where 9 S is the zenith angle for the downwelling radiation impinging on the surface, and 0, is the 
incidence angle for the upwelling radiation leaving the surface in the direction of the 
observation. The integral is over the 2n steradians of the upper hemisphere. The term a 0 , c and 
<jo, x are the co-pol (i.e. polarization equal P) and cross-pol (i.e. polarization orthogonal to P) 
bistatic normalized cross sections, which are functions of unit propagation vectors k, and k s that 
denote the direction of the upwelling and downwelling radiation, respectively. The angle tp s is 
the azimuth of k s relative to kj. These cross sections specify what fraction of power coming from 
k s is scattered into kj. T B d is the brightness temperature of the downwelling atmospheric 
radiation, and T Bcos is the cosmic background plus distant galaxies radiation and equals 3 K (see 
Section 2.2.2). 

The surface reflectivity is given by an equation similar to (32): 

p n/2 2 n 

R p = j sin^J^[cr oc ( k s , k i) + o- ox (k s ,k i )] (33) 

0 0 

The surface emissivity Ep is given by Kirchhoff s law to be 

E p = l-R p (34) 

To numerically compute Tbe from the above theoretical equations, a number of empirical 
models are used as is described in the next sections. 


2.3.2 Empirical Approach to Specifying Atmospheric Terms 

The atmospheric terms are the upwelling and downwelling brightness temperatures Tbu and 
T no. and the atmospheric transmittance r. These are computed by numerically integrating 
equations (28) through (31). For operational processing, we use NCEP profiles of pressure, 
temperature, humidity, and liquid water to find profiles of a(s). The NCEP profiles are 
interpolated to the exact time and location of the Aquarius observations. For the oxygen 
absorption component of a(s) we use the expressions given by Liebe et al. [1992]. For the water 
vapor absorption component of a(s) we use the expressions given by Rosenkranz [1998]. We 
make small modifications to both the Liebe and Rozenkranz formulations based on many years 
of analyzing satellite microwave observations. At 1.4 GHz, the only modification of any 
importance is the non-resonance continuum temperature coefficient which was changed from 
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Liebe’s exponent value of 0.8 to 1.5. This change was based on an analysis of WindSat and 
AMSR-E 7 GHz brightness temperatures. A paper on this is planned for 201 1. 

2.3.3 Empirical Approach to Specifying Emissivity 

The theoretical modeling of the bistatic normalized cross sections is not sufficiently 
advanced to provide sufficiently accurate estimates of Rp and Ep. Instead, we use the 
conventional approach of expressing the emissivity Ep in terms of a specular component (i.e. flat 
surface) and a wind-induced component: 

E p =E 0P + AE p (W,(p r ) (35) 


where relative wind direction <p r is the direction of the wind relative to the azimuth angle of the 
Aquarius look direction. The specular component Eop is computed using the dielectric constant 
model of Meissner and Wentz [2004]. For the wind-induced component, we use an empirical 
model derived from L-band aircraft measurements [Yueh et al., 2010]. Yueh expresses his 
results in terms of an increase in brightness temperature with wind speed: 

AT bv (W) = (0.275-0.0024153 0 +1.4026 -KT 4 Of -2.3326 -10~ 6 0f)w 

{ (36) 

AT bh (IT) = (0.275 + 0.0030010 9 i - 2.5181- 10~ 6 9; -6.9763 -10' 7 0, 3 JW 


where 9, is the Earth incidence angle (deg) and W is the 10-m wind speed (m/s). We convert this 
to a change in emissivity as follows: 


A E P {W,(/> r ) 


A T bp (W) E qp (T s ,S ) 
276.16 E op (276.16,35) 


+ o(w,f) 


(37) 


where the first term is the isotropic wind-induced emissivity and the second term accounts for 
the wind direction variation. Yueh’s measurements were done in cold water around 3°C, and 
hence to convert to emissivity one must divide by 276.16 K. Furthermore, to convert to the 
actual water temperature and salinity of the Aquarius observation, we multiply by the ratio of the 
actual specular emissivity to the specular emissivity corresponding to Yueh’s measurements. 
We assume a nominal salinity of 35 psu for Yueh’s measurements, although this choice has very 
little effect on the calculation. 


The wind direction term for A Ep is small and is yet to be specified. As a placeholder, we are 
assuming the following form for the wind-direction dependence. 

®(W,</) r ) = [(p lP + p 2p 0,.)cos<p,. + (p 3P + p 4 p 9 j ) cos 2(p r ] W (38) 


where the p coefficients are currently set to zero. When the Aquarius Science Team has come to 
agreement on the specification of this term, the p coefficients will be updated. 

For operational processing, the surface temperature Ts is currently the NCEP daily, 0.25° 
product and the wind speed W and direction <p r are the NCEP 6-hour, 1° product. 


2.3.4 Empirical Approach to Specifying Reflected and Scattered Radiation 

The remaining term that needs to be specified is the scattering integral for the downwelling 
radiation: Tppn . A primary reason for computing the space radiation terms separately, as 
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described in Section 2.2, is to facilitate the computation of Tbpq.. In the absence of sharply 
varying space terms, the integral now just contains Tbd+tTb CO s, a small term that varies very 
slowly over angle. The bistatic cross section has a sharp peak in the specular direction at which 
0j = 6 S and tp s = 1 80°, and over this sharp peak Tbd+tT Bc os varies little and can be removed from 
the integral. The remaining integral is simply the surface reflectivity, and hence one has 

TBPn = (T B D + rT Bc J R P (39) 

where T B d and r conveniently have the same slant path angle as T B u- 


2.3.5 Empirical Approach to Handling Land and Sea-Ice Observations 

Although our primary focus is observations over the ocean, the formulation we have 
presented applies equally well to land and sea-ice observations. The only exception is that the 
sea-surface emissivity Ep given by (37) needs to be replaced by a land or sea-ice emissivity, with 
Kirchhoff s law then giving Rp. For observation of mixed surfaces, the following expressions 
are used 


R P = f 


D 

water v Pwater 


fland R Pland + fice R . 


ice Pice 


(40) 


E P T s - ( 1 _ R Pwater ) T water + ( 1 “ R PUmd ) T kind + ( 1 “ R Pice ) T ice ( 4 1 ) 

where / denotes fractional area and T denotes surface temperature. Equations (40) and (41) are 
then substituted into (27). The fractional land f an d is computed from the number of 1/12° land 
pixels that are within an observation footprint and this calculation is based on a static land mask. 
The fractional sea-ice fi ce comes from the NCEP operational sea-ice dataset. The fractional 
water is then f water =1- fi an d ~ fice • The reflectivities for land and sea ice are discussed in 
Appendix A. 


3 Salinity Retrieval Algorithm 

The salinity retrieval algorithm consists of a number of steps that are intended to remove the 
unwanted sources of radiation (galaxy, sun, moon, and Earth’s atmosphere) in order to obtain 
just the Earth’s surface emission term. A simple regression is then use to estimate salinity from 
the surface emission. These steps include: 

1. The radiation received from the galaxy, sun, and moon is subtracted from the Ta 
measurements. 

2. An antenna pattern correction (APC) is done that removes cross-polarization contamination 
and corrects for sidelobes. The APC converts the top-of-the-ionosphere (TOI) antenna 
temperature to a brightness temperature. 

rd 

3. Faraday rotation is removed using the 3 Stokes Tb measurement. This provides the top-of- 
the- atmosphere (TOA) brightness temperature. 

4. The atmospheric components of the TOA Tb are computed using NCEP operational fields 
and are then removed from the TOA Tb. This provides the surface emission by itself. 

5. The salinity value is retrieved from the surface emission using a simple least-squares 
regression. 
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Each of these steps is now described. 

3.1 Removal of the Radiation Coming from Space 

The first step in salinity retrieval is to remove the contribution of space radiation from the 
antenna temperature measurement T\ <mea - Section 2.2 lists the 6 components of space radiation 
and describes how they are computed during operational processing. This subtraction is 
represented by 

T = T -T -T -T -T -T -T (421 

A, earth _dir A ,mea A ,, gal _ dir A,gal_ref A,sun_dir A,sun_ref A, sun _bak A,mon_ref k ' 

where we are now just left with the radiation coming directly from the Earth. 


3.2 Definition of Reference Brightness Temperature 

At this point in describing the retrieval algorithm, we need to define the reference brightness 
temperature that we are seeking to retrieve from T x, ea nh_dir- One possible choice is the Earth 
TOA Tb at boresight. However, we think a better choice is to define the reference brightness 
temperature as an average over the 3-dB Aquarius footprint since this is closer to what is actually 
measured. To be precise, the reference TOA Tb is defined as 




3 dB footprt 


l BE ,toa 


J dA 


3 dB footprt 


(43) 


where the integral is over the 3-dB footprint and the overbar denotes the average over the 
footprint. The 3-db footprint is defined as that part of the Earth integral for which the angle 
between the boresight vector bo and the look vector b is less than 3.07°, 3.17°, and 3.24° for the 
inner, middle, and outer horns, respectively. We further specify that when doing the footprint 
averaging, the incidence angle is held to a constant value 9 i , which is the gain-weighted average 
incidence angle given by 

0, = bought (44) 


where boresight is the boresight incidence angle, and % is 1.00177, 1.00186, 1.00148 for the inner, 
middle, and outer horns, respectively. This adjustment is so small that it probably is not needed 
and one could simply use the boresight incidence angle for 0 j . Thus by definition T I!E um is only 
a function of the Earth scene being viewed, with all antenna pattern effects removed. 

The top-of-the-ionosphere (TOI) reference brightness temperature is denoted by T BEtoi and 
is given by applying Faraday rotation: 

T,e,»,=1'0/)V„ (45) 
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3.3 Estimation of Top-of-the-Ionosphere Brightness Temperature 

Now that we have defined our reference Tb, we must estimate its value given the antenna 
temperature Ta, earth_dir ■ In principle we are seeking the solution of the integral equation given by 
(26): given Ta , ear th_dir what is T BEtoi ? It is not possible to directly invert the integral, and other 
means are required. We use the following method. 

Forward simulations are performed for which a year of simulated Ta , eart h_dir are generated. 
This represents about 25 million observations for each horn at intervals of 1.44 seconds. The 
following regression is then assumed. 




= A T 


A , earth _ dir 


(46) 


where A is a 3 by 3 matrix and is called the antenna pattern correction (APC) matrix and the 
‘hat’ on Tbe, for denotes that this is an estimated quantity as compared to the true value. In 
addition to providing Ta ,earth_dir, the simulation also provides the true value of T BE(oi . The 
standard least-squares method is used to find the APC matrix A that minimizes the variance of 


When deriving the APC matrix A, we excluded simulated observations that contain sea ice 
or that are a mixture of land and water. That is to say, regression (46) is only done for 
observations that are totally land or totally water, i.e., uniform scenes. We did not think it was a 
good idea to train the regression with extremely non-uniform scenes, which pose many 
difficulties when trying to relate antenna temperatures to brightness temperatures. 

Considering all complexities of the integral equation (26), it is not obvious if a simple 
regression l ik e (46) will perform adequately. To evaluate its performance, we compute the rms 

value of T be toi - T BE toi averaged over all the simulated observations. These statistics were done 

separately for all-ocean and all-land observations, but the same APC matrix was used for both. 
Over the ocean, the rms v-pol and h-pol values ranged from 0.05 to 0.08 K. Over land, the rms 
values were considerably higher, ranging from 1 to 2 K. The higher land values are due to the 
land Tb having more variability. A small number of land observations having very large 
variability show errors of 10-20 K, which have a great influence on the rms statistics. These 
results show that the APC regression (46) works well for uniform scenes, l ik e the ocean. Non- 
uniform scenes present a problem for the APC, which is probably unavoidable. 

We experimented with deriving the APC just with ocean observations since this is our 
primary focus. There was a slight improvement when the ocean-only APC was applied to the 
ocean observations, but the improvement was very small compared to the overall 0.05 to 0.08 K 
performance, and we elected to stay with the ocean+land APC, which is trained over the full 
range of T A . 

rd 

The largest error in the APC is with respect to the 3 Stokes for the outer horn, for which 

rd 

the ocean-only observations exhibit an rms of 0.2K. However, the 3 Stokes is just used for 
Faraday rotation correction, and the mapping of this error into the top-of-the-atmosphere (TOA) 
T B is very small. 

The 0.05 to 0.08 K performance of the APC is of some concern, given the demanding 
requirements for Aquarius brightness temperature accuracy. We did some experiments in an 
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attempt to improve its performance, such as using an ocean-only APC, using separate APC for 
ascending and descending orbit segments, and adding higher order terms to (46). These attempts 
helped slightly, but not enough to warrant their implementation. However, with more time and 
analysis, we do think that there is the potential for improving the APC performance. 

3.4 Correction for Faraday Rotation 

The next step in the retrieval is a correction for Faraday rotation. Yueh [2000] gave a 
detailed description of the Faraday rotation problem and proposed measuring the 3 ld Stokes 
parameter to remove its effect, as we will now discuss. Faraday rotation is removed by applying 
the following expression. 

T K ,„='l'H)t BI ,„, (47) 

If one assumes that the 3 ld and 4 th TOA Stokes from the Earth are negligibly small, one may set 
the T 1{E toa 3 ld Stokes to zero and obtain the Faraday rotation angle: 




1 

= — arctan 
2 


( 


V 


T 

' BE.Ioi , 3 

f 

± BE,toi,2 


\ 


J 


(48) 


where the last subscript denotes Stokes number. Substituting (48) into (47) one obtains the TOA 
T b 


f =f 

± BE,toa,\ x BE,toi,\ 


1 BE,toa,2 


* 

2 

r a, n 

T 

1 BE,toi,2 1 

+ 

T 

1 BE,toi ,3 


“|2 


(49) 


Wentz [2006] analyzed the error between the estimated T BE toa and the true T BE toa . The rms 

difference was about 0.04 K for ocean v-pol and h-pol observations. The Wentz [2006] results 
used an antenna gain pattern based on a JPL computer model. Since then, the computer model 
antenna pattern has been replaced by the scale-model antenna pattern, which has considerably 
larger cross-polarization terms. These large cross -polarization terms result in larger error when 

doing the APC. The T BE tou - T BEtoa rms value is now 0.05 to 0.08 K, nearly all of which is due to 
the APC error as mentioned in the previous section. 


3.5 Removal of Atmospheric Effects 

In this subsection we switch over to using the modified Stokes parameters (v-pol and h-pol) 
rather than the classical Stokes parameters for T BEtofl , and these two components are denoted 

simply as T B e,p, where subscript P denotes polarization (either V or H). This step in the retrieval 
algorithm is intended to remove the effect of the atmosphere. This is done by inverting (27) to 
yield surface emission T B E,sur ■ 


T -T 

= E T — BE - P BU 


s 


-T, 


bp n 


(50) 


Substituting (39) in (50) gives 
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BE. P. sur P S 


= E.T, = — ~(T bd + tT b )R, 


Noting that Rp =1 - E P , (51) can be rewritten as 


T -T 

± BE.P ± BU ( rp rp \ 

(T bd + tT B c J 


T S + ( T BD + tT Bcos ) 


(51) 


(52) 


where the bracketed term is the surface emissivity. 

The atmospheric terms for the upwelling and downwelling brightness temperature Tbu and 
T no, and the atmospheric transmittance r are computed by numerically integrating equations (28) 
through (31). For operational processing, we use NCEP profiles of pressure, temperature, 
humidity, and liquid water to find profiles of a(s). The NCEP profiles are interpolated to the 
exact time and location of the Aquarius observations. See Section 2.3.2 for more details. The 
surface temperature Ts currently comes from the NCEP daily, 0.25° product. After launch, other 
sources for Ts will be considered such as AMSR-E and WindSat microwave-only retrieval of sea 
surface temperature. T Bc0 s is a constant 3 K. Given these ancillary data, (52) can be calculated to 
yield the surface emission brightness temperature. Note that this calculation of Tbe.p.suv is done 
for all surface types (ocean, land, sea-ice, and mixed). 


3.6 Estimation of Salinity Given Surface Emission Brightness Temperature 

Over the ocean, the surface emission TpE.p.sur depends on surface temperature T s , surface 
salinity S, surface roughness, and Earth incidence angle For now, we are parameterizing the 
surface roughness effect simply as a function of the wind measured 10 meters above the surface, 
where W denotes the speed (m/s) and (p r denotes the relative wind direction. The wind direction 
effect given by (38) is first removed from the surface Tb. 

= (53) 

where the prime sign denotes that the wind direction effect is removed. 

Wentz and Yueh [2004] described the algorithm for estimating sea-surface salinity from the 
surface emission T' BE p sur . The algorithm takes the form 


S = S 0 ( 0 i d S )+S l {0„t S ) TpEy^ur + S 2 {Opts) T BE,H .sur + S 3 {0„t s )W 


(54) 


where 0, is the incidence angle, ts is the sea-surface temperature in °C ( T s = ts + 273.15), and W 
is sea-surface wind speed (m/s). The s coefficients are functions of Q, and ts and are in tabular 
form. The algorithm is trained by deriving a set of s coefficients that minimize variance between 
the salinity S given by (54) and the true salinity. Separate sets of 5 coefficients are found for 25 1 
incidence angles going from 25° to 50° in 0.1° steps and 451 ts values going from -5°C to 40°C 
in 0.1 °C steps. These 251x451 tables are interpolated to the specified value for 6, and ts. 
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The algorithm training requires that T' HE p sur be computed over the full range of S. T s , W and 

0i using equations (35) through (38). Note that this is the only point in the retrieval algorithm 
where the specification of the dielectric constant of sea water [ Meissner and Wentz, 2004], and 
the wind-induced emissivity [Yueh et al., 2010] comes into play. Using a different dielectric 
model or a different wind model would result in a different set of s coefficients. 

For a given Q, and ts, the relationship between surface brightness temperature and salinity is 
close to linear, but not quite. If necessary, we will supplement equation (54) by adding a second 
stage algorithm that will remove any error caused by the deviation from linearity. This type of 
two-stage algorithm has proven very effective in handling the non-linearity characteristic in the 
AMSR-E retrievals. However, any such fine tuning should wait until the L-band wind model 
(i.e., equations 36-38) is known more precisely. We expect that about 6-months after launch an 
improved wind model will be developed, and we can then take a closer look equation (54). 

Note that in training the algorithm, a different amount of noise can be applied to the v-pol 
and h-pol Tb to account for the actual error in these Tb retrievals. Also, an analogous algorithm 
can be developed that uses the sum of v-pol and h-pol Tb to eliminate residual Faraday rotation 
error. The salinity retrievals from the V+H algorithm can be compared to those coming from the 
dual-polarization algorithm to test if there is any residual Faraday rotation error. This type of 
experimenting will be done after launch. 

Another alternative to (54) is to subtract the wind-induced Tb correction A Tbp(W) from 
Tbe.p.sut prior to the salinity retrieval. This is the approach used for the algorithm derived during 
the WISE experiments [Camps et al., 2002]. However as said above, it is probably best to wait 
until we have the improved post-launch wind model before trying these different alternatives. 

The salinity retrieval is very sensitive to errors in Ts and W. The pre-launch version of the 
algorithm currently uses the NCEP daily, 0.25° product for Ts and the NCEP 6-hour, 1° product 
for W. After launch, other sources for Ts and W will be considered. Microwave retrievals of sea 
surface temperature from AMSR-E and WindSat could possibly prove to be more accurate than 
the NCEP sea surface temperatures. Also, at some point, we expect to start using the wind speed 
derived from the Aquarius microwave scatterometer. 


3.7 Correction for Land Entering the Antenna Sidelobes 

One option that is being evaluated is the removal of sidelobe contamination when the 
Aquarius observation is close to land. Simulations have shown that the salinity retrieval 
degrades quickly as the footprint approaches land (closer than 400 km). This land-contamination 
error, which is shown in the Wentz [2006] study, is defined as 

= T B E,toa ~ T BE ,toa (5^) 

which is the difference of the estimated TOA Tb minus the true TOA Tb coming from the 
simulation. We plan to use the Aquarius simulator to produce tables of AT BE(oa . These tables 
will be stratified according to latitude (0.25° increment), longitude (0.25° increment), 
ascending/descending (1 or 2), month (12), horn (3), and polarization (2). The monthly variable 
is required to account for seasonal variability in the land brightness temperatures and the 
variability of sea ice. (Sea ice will be particularly challenging due to interannual variability.) 
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The land-correction tables will have the form Tbjand_correction(1440,720,2,12,3,2). In principle 
these table values could be subtracted from the estimated TOA Tb in order to remove the land 
contamination. Possibly, this might provide usable salinity retrievals closer to land: maybe to 
within 200 km of land. 

Other methods for the removal on land contamination are also being investigated. One 
method for sidelobe removal that has been applied to radiometers on JASON and ENVISAT is to 
divide the integration (e.g. in equation 26) into additional sectors such as the on-Earth sidelobes 
to obtain a mainbeam brightness temperature, and to eventually deconvolve the mainbeam 
brightness temperature [Brown, 2005]. 


4 Remaining Work before Launch 

The biggest challenge we face before launch is improving our modeling of the reflected 
galactic radiation (Section 2.2.2). Currently, we are only accounting for galactic reflections in 
the specular direction. In truth, bistatic scattering from a rough ocean will scatter galactic 
radiation from many different directions into the mainlobe of the antenna. Between now and 
launch, we plan to include bistatic scattering in our simulator. This is a formidable undertaking 
because it means putting a double integral into the kernel of the Ta integral equation (4). 
However, we have done some initial testing, and the problem is tractable although the processing 
time is considerable. We estimate it will take 32 processors about 2 months to do the necessary 
calculations. No changes to the salinity retrieval algorithm will be required other than updating 
the galactic reflection tables tagal_ref_tab(1441 ,1441 ,3,3,5). These tables have already been 
parameterized in term of wind speed in anticipation of this update. 

Another issue that needs revolving before launch is the specification of the wind directional 
term for O ( W . </> ) given by (38). Again this should not require changes to the retrieval 
algorithm other than specifying the p coefficients, which are currently set to zero. 

Finally, we want to implement the land correction tables discussed in Section 3.7. Whether 
this is done before or after launch needs to be decided. This change will require a small addition 
to the current code. 

Finally, research is underway at the George Washington University and at the Polytechnic 
University of Catalonia to refine the Klein and Swift [1977] model function for emissivity at L- 
band that relates the brightness temperature to sea surface temperature and salinity. After 
launch, the salinity retrievals using the Meissner-Wentz dielectric constant model can be 
compared with those using the alternative dielectric constant model. This type of experimenting 
is done by recomputing the s coefficients in (54) using a different dielectric constant model. 
Similar experimenting can be done with different wind-induced emissivity models. 
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6 Appendix A. Land and Sea-Ice Reflectivities 

This appendix will be done soon. It will contain expressions for computing the reflectivity of 
land and sea ice. 
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